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ABSTRACT 

Future large ensembles of time delay lenses have the potential to provide interesting cosmological 
constraints complementary to those of other methods. In a flat universe with constant w including 
a Planck prior, LSST time delay measurements for ~ 4, 000 lenses should constrain the local Hubble 
constant h to ~ 0.007 (~ 1%), tide to ~ 0.005, and w to ~ 0.026 (all 1-a precisions). Similar constraints 
could be obtained by a dedicated gravitational lens observatory (OMECA) which would obtain precise 
time delay and mass model measurements for ~ 100 well-studied lenses. We compare these constraints 
(as well as those for a more general cosmology) to the "optimistic Stage IV" constraints expected from 
weak lensing, supernovae, baryon acoustic oscillations, and cluster counts, as calculated by the Dark 
Energy Task Force. Time delays yield a modest constraint on a time- varying w(z), with the best 
constraint on w(z) at the "pivot redshift" of z ~ 0.31. Our Fisher matrix calculation is provided 
to allow time delay constraints to be easily compared to and combined with constraints from other 
experiments. We also show how cosmological constraining power varies as a function of numbers of 
lenses, lens model uncertainty, time delay precision, redshift precision, and the ratio of four-image to 
two-image lenses. 

Subject headings: cosmological parameters - dark matter — distance scale — galaxies: halos - 
gravitational lensing — quasars: general 



1. INTRODUCTION 

The HST Key Project relied on 40 Cepheids to con- 
strain Hubble's constant H to 11% (jFreedman et alJ 
2001). The first convincing measurements of the ac- 
celerating expansion rate of th e universe ( sugge sting 
the existence of dark en ergy) by iRiess et al.l (|1998f ) and 
iPerlmutter et al.1 (|1999f ) required 50 and 60 supernovae, 
respectively. So far, time delays have only been reliably 
measured for ~ 16 gravitational lenses, thanks to dedi- 
cated le ns monitoring from cam paigns such as COSMO- 
GRAIL (jEigenbrod et al.ll2005h . Yet recent analyses of 
10-16 time delay lenses already claim to match or surpass 
the Key Project's 11% p recision on H (ISaha et al.ll2006l : 
lOgurill2007l : lColesj |2008V Future surveys promise to yield 
hundreds or even thousands of lenses with well-measured 
time delays, which will enable us to obtain much tighter 
constraints on Hq as well as constraints on other cosmo- 
logical parameters. 

To date, most efforts have focused on studies of indi- 
vidual time delay lenses. In theory, one might be able 
to control all systematics and constrain Hq unambigu- 
ously given a single "golden lens" . Such a lens would 
have a sufficiently simple and well-measured geometry. 
The closest to a g olden lens may be B1608+656. In 
iSuvu et ail (|2009b), the authors claim all systematics 
have been controlled to 5%. A new estimate for H based 
on this lens is forthcoming (jSuvu eTaD[2QQ93). 

Historically, analyses of individual lenses hav e yielded 
varyi ng answers for Hq (see the Appendix of UacksorJ 
120071 for a recent review). This can be attributed to 
two factors, both of which, it appears, are now being 
overcome. 

The first factor is simple intrinsic variation in lens 
properties (especially mass slope) and environment (lens- 
ing contributions from neighboring galaxies). Consider 



the following estimate from a simple empirical argument. 
If statistical uncertainties on Hq decrease as 1/ y/~N (as- 
suming systematics can be controlled), and the current 
uncertainty from 16 lenses is ~ 10%, then the uncer- 
tainty on a single lens might be ~ 40%. Thus, assuming 
h = 0.7 (where H = lOO/i km s^ 1 Mpc -1 ), individual 
lenses may be expected to yield a wide range of h = 0.42 
- 0.98 (1-<t) . (We will revisit these assumptions in this 
work.) 

The second factor in the wide range of reported Ho 
values is that different analyses have assumed different 
mass profiles to model the lenses, including isothermal, 
de Vaucouleurs, and mass follows light. There is sub- 
stantial weight of evidence that galaxy lenses are roughly 
isothermal on aver age, at least within app roximately the 
scale radius (e.g.. iKoopmans et al.l |200"6T) . Theoretical 
work supports this idea, showing that a wide range of 
plausible luminous plus dark matter profiles all combine 
to yield roughly an isothermal profile at the Einstein ra- 
dius, though the slope may deviate from isothermal be- 
yond that radius ijvan de Ven et al] l2009). 

In recent years we have witnessed a steady in- 
crease in the number of strong lenses disco vered by 
searches such as C LASS (IMvers et al l [2001. SLACS 
(iBolton et al.ll2006h . SL2S ICabanac et al.ll2007lL SOL S 
(|lnada et all 12008). HAGG LeS dMarshall et"afll2009bft . 
and search es of AEGIS dMqustakas et al.l l2007f) and 
COSMOS (|Faure et all MMf. Based on this ex- 
perience, we can expect that future s urveys such 
as Pan-STARRS 1 (IKaiserl 12001 . LSST 2 (jlvezic et al.1 

1 The Panoramic Survey Telescope & Rapid Response System, 
http: / /pan-starrs. ifa.hawaii.edu 

2 The Large Synoptic Survey Telescope, http://www.lsst.org 
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l2008h . JDEM / IDECS 3 , and SKA 4 l|Laziol f2008h will 
yield an expl osion in the number o f strong lenses 
known (e.g., IKoopmans et al.l 120041 : IFassnach t et al.l 
120041 Marshall et all 120051 ). Prospects for using these 
lenses to constrain the nature of dark matter over 

the course of the nex t decade were present ed in 

Moustakas et alj (120091) . IKoopmans et all (|2009al) . and 
Marshall et alj (|2009al) . 



It is reasonable to expect that time delays will be 
reliably measured for large numbers of these lenses, 
whether through repeated observations in surveys (Pan- 
STARRS and LSST), auxiliary monitoring, and/or 
through tailored spec ific missions such as OMEGA 
(Moustakas et al. 2008). Increased sample size, improved 
lens model constraints, and higher precision redshifts and 
time delay measurements will all improve constraints on 
Hq and other cosmological parameters, as we present be- 
low. 

A more precise measurement of Hq will yield tighter 
constraints on both the dark energy equation of state 
parameter (w) and the flatness of our universe (flk), 
indep endently of the results of fu t ure dark energy sur- 
veys ([Blake et al.l 120041: iHul 120051 : lAlbrecht et"alJl2006t 
lOllinel 120071 ). To this end, the SHOES Program (Su- 
pernovae and Hq for the Equation of State) has ob- 
tained new observations of supernovae and C epheid vari- 
ables with reduced systematics. Recently, [Riess et alj 
(2009) published a redetermination of Hq = 74.2 ± 
3.6km s _1 Mpc -1 , or 5% uncertainty including both 
statistical and systematic errors. Their Hq determi- 
nation plus WMAP 5-year data alone constrain w = 
— 1.12 ± 0.12 ( assum ing constant w). 

iRiess et~a l. (2009) also make the following important 
point that bears repeating. The seemingly tight con- 
straints on Hq derived from CMB + BAO + SN exper- 
iments are in fact predictions or inferences of Hq given 
those data and a cosmological model. They are no substi- 
tute for direct measurement of Hq such as that presented 
i n their work o r the HST Key Project. 

Oiling ( 2007) reviews several methods with the poten- 
tial to directly constrai n Hq. Water masers, for exam - 
ple, hold much promise (|Braatz et alj2008tlBraat92009| ). 
Time delays and water masers both yield direct geomet- 
ric measurements of the universe to the redshifts of the 
observed sources (z ~ 2 or greater for time delay lenses), 
bypassing all distance ladders. 

Time delays do not simply constrain Hq. To first order, 
each time delay is proportional to the angular diameter 
distance to the lensed object and thus inversely propor- 
tional to Hq. An additional factor involves a ratio of two 
other distances - from observer to lens and from lens 
to source. All three of these distances have a complex 
(though weaker) dependence on the other cosmological 
parameters (fl m , fide, flk, wo, w a ) which contribute to the 
expansion history of the universe. 

Most time delay analyses ignore this weaker depen- 
dence on (f2 m , fide, Ofe, wo, w ), in effect assuming these 
parameters are known perfectly. In this paper we show 
how relaxing this "perfect prior" increases the uncertain- 
ties on H . As dark energy surveys endeavor to place con- 
straints on w and the flatness of our universe fl^, we must 

3 The Joint Dark Energy Mission, http://jdem.gsfc.nasa.gov 

4 The Square Kilometer Array, http://www.skatelescope.org 



study how time delays can contribute to these constraints 
without assuming the very parameters we would like to 
constrain. In this work we also study the ability of large 
time delay ensembles to constrain (f2 m , fide, flk, wq, w a ). 

The idea to use time de l ay len ses to measure Hq was 
first proposed by iRefsdall ([1964) . Strong gravitational 
lenses are elegant geometric consequences of how light 
travels through the universe while grazing massive galax- 
ies. When the line of sight alignment is very close, light 
takes multiple paths around the curved space of the lens. 
These paths form multiple images, and the light takes 
a different amount of time to travel each path. Light 
passing closer to the lens is deflected by a larger angle 
(increasing its path length) and experiences a greater rel- 
ativistic time dilation, further delaying its arrival. If the 
source flares up, or otherwise varies in intensity (e.g., 
if it is an active galactic nucleus, or AGN), we can ob- 
serve these "time delays" between or among the images. 
These time delays are functions of the angular diameter 
distances between the source, lens, and observer, as well 
as the properties of the lens itself. 

The ability of time delays to constrain o ther cosmolog- 
ical p arameters has also been explored. iLewis fe Ibatal 
(2002) explored various combinations of (fl m ,fl\) in a 
flat universe and various (wo,w a ) for fixed (fl m ,fld e )- 
Most notably, they calculated constraints on (h, w) from 
ensembles of lenses assuming constant w and (fl m ,fl\) 
= (0.3, 0.7), finding that h and w would not be strongly 
constrained. We show that the addition of a Planck 
prior improves these constraints considerably. iLinderl 
(2004) investigated constraints on the dark energy pa- 
rameters (wq,w„) from various methods, touting the 
complimentarity of strong lensing to that of other meth- 
ods. However, they concede that the unique positive 
correlation in strong lensing (wo,w ) constraints evapo- 
rates whe n including degeneraci e s othe r cos mological pa- 
rameters. iMortsell fc Sunessonl (|2006f) and iDobke et all 
(2009) examined the constraints that large ensembles 
of lenses might place on Hq and fl\ = 1 — fl m (as- 
suming a flat universe). Below we present the first 
full treatment of the cosmological constraints expected 
on (h, fl m , fide, flk, wo, w a ) from ensembles of time delay 
lenses including various priors. 

Lens statistics from well-controlled searches for 
strongly-lensed s ources have also been used to constrain 
cosmology (e.g.. IChael 120071 : lOguri et all 12001 . If time 
delays can be obtained for the lenses in such a sample, 
the lens statistics and time delays might combine to yield 
tighter cosmological constraints. This potential is not ex- 
plored in this work. 

Cosmological constraints can also be obtained from 
symmetric strong lenses for whi ch velocity disper- 
sions have been measured (e.g., IPaczvnski fc Gorskil 



19811 fF utamas e fc Hamava!ll999l4Yamamoto et alJl200H 
Lee fc Nd I2007D . Assuming an isothermal model, the 
measured velocity dispersion determines the Einstein ra- 
dius solely as a function of cosmolog y (given redshifts 
measu red to the lens and source). lYamamoto et al.l 
(|2001| ) studied the future potential for this method to 
constrain cosmology using a Fisher matrix analysis. 

The reader is invited to skip ahead to our results in 
$51 where cosmological constraints expected from time 
delays (according to our calculations) are compared to 
those expected from other methods (weak lensing, super- 
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novae, baryon acoustic oscillations, and cluster counts). 
Table [2] summarizes the assumed priors including a guide 
to specific sections and figures. 

The remainder of our paper is organized as follows. In 
$2] we provide the time delay equations and discuss how 
cosmology is derived from observed time delays. We de- 
fine the quantity T c (h, Sl m , Side, Sl k , w , w a ; z L , z s ) which 
time delays are capable of constraining. In <|3]we estimate 
the constraints on Tc expected from future experiments. 
(A more detailed analysis of l ensing simulations is pre- 
sented in a companion paper iCoe &: Moustak as 2009a, 
hereafter Paper I.) In S|4]we illustrate the dependence of 
Tc on cosmological parameters (h, fi m , Side-, wo, w a ). In 
as highlighted above, we give projections for time 
delay constraints on (h, Side, Slk, wo, w a ) and compare to 
other methods. Systematic biases are discussed in $6] and 
their impact on our ability to constrai n cosmology is an- 
alyzed in another companion paper (jCoe &: Moustakaa 
2009c, hereafter Paper III). Finally we present our con- 
clusions in $7l 

We assume all constraints to be centered on the concor- 
dance cosmology h = 0.7, Sl m — 0.3, fide = 0.7, Slk = 0, 
wq = —I, and w a = 0, where Hq = lOO/i km s _1 Mpc -1 . 

2. COSMOLOGICAL CONSTRAINTS FROM TIME DELAYS 

2.1. Time Delay Equations 

A galaxy at redshift zl strongly lenses a background 
galaxy at redshift zs to produce multiple images. Either 
two or four images are typically produced. 5 We refer 
to these cases as "doubles" and "quads", respectively. 
The lensing effect delays each image in reaching our 
telescope by a different amount of time, given by 



\e-(3\ 2 - 



(i) 



(e.g., Blandfo rd fc Naravanl 119861 ) with terms defined 
below. The factors in the time delay equation can be 
grouped into a product of two terms: 



The first factor, 



At = T c T c . 



t c = ii±^p, 



(2) 
(3) 



is a function of cosmology and the lens and source 
redshifts, Zi and zs- The second factor, 



Tr 



(4) 



is a ratio of the angular-diameter distances from ob- 
server to lens Dl = Da(0,Zl), observer to source 
Ds = Da(0,zs), an d lens to source D L s = Da{z l ,z s ). 
Angular-diameter dis tances are calculated as follows 
dFukugita et al.1 119921 filled beam approximation; see 
aiso[Hogg 1991: 



D A (zi,z 2 ) 



E A = 



c E A (z\,z 2 ) 



Hn 1 



Z2 



\/\£lk\E\ 



(6) 



(7) 



where sinn(w) = sin(it), u, or sinh(w) for an open, flat, 
or closed universe respectively (fife < 0, Slk — 0, or 
Qk > 0)- The curvature is given by fife = 1 — (fi m + ^a)j 
while 



dz' 



(8) 



The normalized Hubble parameter E(z) can have dif- 
ferent expressions depending on the cosmology assumed: 



H{z) 



(9) 



= y/n m (i + z) 3 + n*(i + z) 2 + 

= ^r» m (l + Z)3 + n fc (l + Zf + Oie(l + Z)3(l+w) 

= ]j--- + + z) 3 ( 1+w »+ w «> exp (^^)- 

Here we have progressed from a universe with a cosmo- 
logical constant J7a to one with dark energy with an 
equation of state p = wp. In the last line, the last term 
has been rewritten in terms of an evolving dark energy 
equation of state 



w = w + w Q (l — a) 

z 



■■ w + w a 



1 + js 



(10) 
(11) 



a common parame t rizati on first in t roduc ed by 
IChevallier fc Polarskil (pOOl and ILinderl (|200l . The 
universe scale factor a = (1 + 
We next define the dimensionlcss ratio 



is a function of the projected lens potential <f>, the source 
galaxy's position on the sky /3, and the image positions 

e. 

We concentrate on the cosmological dependence of Tc . 
The factor 



V = 



DlDs 
Dls 



(5) 



5 An additional central demagnified image is also produced by 
every lens with a central mass profile shallower than isothermal. 
Such images are rarely bright enough to be detected, thus we ignore 
them throughout this work. 



E = 



ElEs 



E 



LS 



(12) 



with factors defined similarly to those above for D 4: 
E L = E A (0,z L ), E s = E A (0,z s ), E LS = E A {z L ,z s ). 
We find that many factors cancel, and Tc simplifies to: 



T c = 



£(Q m , Side, flk, Wp, W a ) 

Ho 



(13) 



We see here clearly that time delays (At = TcXc) scale 
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inversely with Hq. There is also a complex though weaker 
dependence on the other cosmological parameters as em- 
bedded in £. 

2.2. Deriving Cosmology from Time Delays 

Given observed time delays At and assuming a lens 
model (and thus Tc), one can obtain measures of Tq. 
These measures will have some scatter due to both ob- 
servational uncertainties and deviations of the lens from 
the assumed model. 

Recent studies suggest that galaxy lenses, on average, 
have roughly isothermal profiles within the Einstein 
radius (see SJTJ). Deviations from this simple description 
include variation in lens s lope, e xterna l shear, mass 
sheets, and substructure. IQguril ((2007) parametrized 
the deviations as the "reduced time delay", the ratio 
of the observed time delay to that expected due to an 
isothermal potential in a given lens: 



A7 



Ar if 



(14) 



In our notation, these observed deviations are due to 
deviations in the lens model: 



Tr 



T 



C.iso 



(15) 



By assuming an isothermal model (Tc = T^iso), these 
deviations get absorbed into the derived cosmology: 



Tr 



Tr 



C.truc 



(16) 



where 7c, true 1S ^ rie true cosmology. For example, a lens 
which is steeper than isothermal yields S£ > 1; thus 
when assuming an isothermal model (S£ = 1), we derive 
Sc > 1 (since EE = EEeEEc). In traditional analyses as- 
suming fixed £, Sc > 1 would simply yield a low h. This 
approximation is adequate for small samples of lenses 
but not for the large samples to come in the near future 

(9533). 

Similarly, observational uncertainties affecting At are 
absorbed into the derived cosmology. In this paper, we 
study how observational and intrinsic (lens model) un- 
certainties combine to yield scatter in the observed At. 
We will assume these measurements yield Tc with the 
correct mean but a simple Gaussian scatter and explore 
how this propagates to Gaussian uncertainties on cosmo- 
logical parameters. 

In practice we do not expect EE£ and measurements 
of At to have Gaussian scatter, but these serve as use- 
ful approximations. The true expected -P(EE) from time 
delay measurements and methods for handling these dis- 
tributions are studied in IQguril (|2007T ) and Paper I. 

3. CONSTRAINTS ON T c FROM FUTURE EXPERIMENTS 

3.1. Extrapolating from Current Empirical Results 

Recent studies have constrained Tc to ~ 10% using 
time delays, where Tc encodes all of the cosmological de- 
pendencies f H2.1[) . Constraints on Tc have generally been 
interpreted to be equivalent to direct constraints on h. 
This assumption is reasonable for current sample sizes, 



but will nee d to be revise d in the future ( H5.4.ip . Us- 
ing 16 lenses. lOgurl (120071) obt ain h = 0.70 ± 0.06 ( stat.) 
Similar studies by iSaha et all (2006) and iColesI (j2008h 
using a different method obtain similar constraints us- 
ing 10 and 11 lenses, respectively. The latter finds 



0.71 



+0.06 
-0.08- 



We will adopt the IQguril (|2007f ) uncertainty of 8.6% 
with 16 lenses as the "current" uncertainty in Tc- 6 We 
note that the time delay uncertainties in this sample are 
roughly and broadly scattered about A(At) = 2 days. 7 

We can improve on these Tc constraints in three ways: 
obtaining larger samples of lenses, better constraining 
our lens models, and obtaining more precise time delay 
measurements. As we explain below, we expect future 
surveys such as Pan-STARRS and LSST to improve on 
the sample size while the lens model and time delay un- 
certainties will remain about the same. These surveys 
will have to contend with a lack of spectroscopic red- 
shifts for most objects, but the gains in sample size will 
more than compensate. Similarly tight constraints on Tc 
could also be obtained by studying relatively fewer lenses 
in great detail, as we discuss below. 

Here we consider statistical uncertainties only, with 
systematics to be discussed in ^SJ We will assume that 
all other things being equal, increasing our sample size 
beats down our errors by y/N for N lenses. This assump- 
tion is borne out well by our detailed simulations (Paper 
I), for the case of no systematic uncertainties. 

Based o n the curren t constraint of STc « 8.6% from 
16 lenses (jOgurill2007h . we project that simply increas- 
ing the sample of lenses would produce constraints of 
ST C ~ 34%/VJV. We will define this as the uncer- 
tainty from lens models and time delay measurements: 
SEct ~ 0.344. Photometric redshifts would degrade 
these constraints as estimated below (<|3~3|). 

3.2. Future Surveys 

Pan-STARRS and LSST will both survey the sky re- 
peatedly, opening the time domain window for astronom- 
ical study over vast solid angles. Pan-STARRS 1 (PS1) 
has recently begun its 3n survey, repeatedly observing 
the entire visible sky to ~23rd magnitude every week 
over a 3-year period. LSST promises similar coverage 
and depth every 3 nights with first light scheduled for 
2014. 

These surveys will reveal many time- variable sources, 
among them gravitationally-lensed quasars. The persis- 
tent monitoring over many years should yield time delays 
"for free" for many strongly-lensed quasars. Simulations 
(M. Oguri 2009, private communication) show that Pan- 
STARRS 1 and LSST are expected to yield - 1, 000 and 
~ 4,000 strongly-lensed quasars with quad fractions of 
19% and 14%, respectively. 

We will assume that these surveys will measure time 
delays to about 2-day precision, or similar to that of 

6 The Oguri (2007) simulations initially suggested an uncertainty 
of ~ 4% in Tc. However jackknife resampling of the data revealed 
the true uncertainty to be twice as much. Under-prescribed shear 
in the simulations was cited as a potential cause for the discrepancy. 

7 We adopt a notation in which "A" refers to uncertainties with 
units and "<5" to fractional uncertainties. Thus a time delay of 
20 days measured to 2-day precision has A (At) = 2 days and 
<5(Ar) = 0.1. 
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redsh ift 



Fig. 1. — Distributions of lens and source redshifts used in this 
paper. These Gaussian distrib utions (z^ = . 5 ± .15, zg = 2.0 ± 
0.75; zg > zl) were used by Dobke ct al. (20Q3) as reasonable 
approximations for near-future missions including LSST. 

our current sample of time delay lenses. This is con- 
sis tent with predic t ions b ased on detailed simulations 
bv lEigenbrod et aT] (|2005| ) which study factors including 
survey cadence, object visibility, and the complicating 
effects of microlensing. We note this estimate may be 
a bit optimistic for PS1 with its slower sampling rate 
compared to LSST. 

The expected redshift distributions of the lenses and 
sources can be roughly approximated by the Gaussian 
distributions zl = 0.5 ± 0.15 and zs = 2.0 ± 0. 7 5 wit h 
z s > z L (Fig. [J), as adopted by iDobke et all ([20091 ) . 
Obviously the two distributions will be correlated, but 
we approximate them as being independent. 

As surveys attain fainter magnitude limits, it is be- 
lieved that the magnification bias enjoyed by quads will 
be diminished. Future surveys are thus expected to yield 
lower quad fractions (~ 19%, 14%) than the current sam- 
ple of time delay lenses (6 / 16 = 37.5%). This might 
improve the expected constraints on 7c from future sur- 
veys as quads have been shown to yield time dela ys with 
more scatter and thus less reliable estimates of 7c (jOguril 
2007, Paper I). 8 However, we find this to be mitigated 
by the fact that quads yield multiple time delay measure- 
ments (one for each pair of images), while doubles only 
yield a single At measurement. Based on our detailed 
simulations and analysis (Paper I), we find quads and 
doubles to have approximately equal power to constrain 
7c • This simplifies our analysis; the quad-to-double ratio 
need not be considered when estimating STq for a given 
experiment. To allay any concern, we stress that this 
assumption actually makes our estimates of 5Tc more 
conservative for future surveys which have lower quad 
fractions than the current sample. 

For each double or quad, image pairs can be further 
classified by their geometry. For example, image pairs 
with small opening angl es are found to yield larger scat- 
ter in At (|Qguril l2007. Paper I). Detailed analyses in 
these papers quantify these scatters, enabling a well- 
informed prior -P(3) to be placed on each image pair as a 
function of geometry. The details are unimportant here 

8 This is believed to be due to the fact that some of the fac- 
tors (especially external shear) which cause scatter in 3 also raise 
the likelihood that a lens will produce quad images rather than a 
double. 



though we hav e made use o f the constraint this analysis 
has put on T c (|Ogurill2007f ). 



3.3. Photometric Redshift Uncertainties 

Currently all lenses which have reliable time delay mea- 
surements also have spectros copic redshif ts measured for 
both lenses and sources (e.g. JOgurill2007l ). The telescope 
time required to obtain spectroscopic redshifts is gener- 
ally a small fraction of that required to obtain accurate 
time delays, so the extra investment is worthwhile. 

Future surveys which repeatedly scan the sky, how- 
ever, will yield time delays for many more lenses than 
may be followed up spectroscopically. For these lenses 
we will have to rely on photometric redshift measure- 
ments. These uncertainties will degrade the constraints 
possible on the cosmological parameters. 

Photometric redshift uncertainties for the lenses (typ- 
ically elliptical galaxies at zl ~ 0.5) are expected to be 
Az L ~ 0.04(1 + z L ), similar to th at found in the CFHT 
Legacy Survey (Ilb ert et al1l2006f) . Redshift uncertain- 
ties for the lensed sources (quasars) are expected to be 
somewhat higher. We will adopt Azs ~ 0.10(1 + zs), 
roughly that f ound in the ana lysis of ~one million SDSS 
quasars ([Richards et al J 120091 ) . 

Obtaining photometric redshifts in ground-based im- 
ages will often be complicated by cross-contamination of 
flux among the lens and multiple images. Yet improved 
photometric redshift techniques are also be ing developed 
with LSST in mind ([Schmidt et al.ll2009f ). so it is per- 
haps too early to say whether our estimated redshift un- 
certainties are too optimistic or pessimistic for a future 
ground-based survey. Some of the most common catas- 
trophic redshift degeneracies can clearly be avoided by 
considering the observed image separations, time delays, 
etc. Most obviously, th e common degen eracy between 
z ~ 0.2 and z ~ 3 (e.g., iCoe et al.1 (2006) can be neatly 
averted since a lens at z ~ 3 or a source at z ~ 0.2 would 
clearly stand out. 

Assuming the above redshift uncertainties, we now de- 
termine how these propagate into uncertainties on 7c. 
For simplicity, let us assume that redshift uncertainties 
are Gaussian. Let us further assume that uncertainty 
in S scales linearly with redshift uncertainty. (This 
is approximately true for reasonable uncertainty levels 
Az<0.2.) 

Using equations [7] - [13l we find for a typical lens- 
source combination with (zl,Zs) = (0.5,2.0), that lens 
and source redshift uncertainties translate to SEz L ~ 
2.75Azl and SEz s ~ — 0.16Azs, respectively. Given the 
above redshift uncertainties, these evaluate to SEz L ~ 
0.16 and SEzg ~ 0.05. These relations are strong func- 
tions of redshift and become catastrophic for sources very 
close to the lens. We plot this behavior in Fig. [2l If ac- 
curate and precise redshifts are not available, we must 
concentrate our analysis on systems with high separa- 
tion in redshift between the lens and source. 

For a lens ensemble with Gaussian redshift distri- 
butions zl = 0.5 ± 0.15 and z$ = 2.0 ± 0.75, we find 
SEz L ~ 0.175 and SEz s ~ 0.028. To calculate these 
uncertainties, we sum the x 2 °f individual lens-source 
combinations, weighting by the probability Pi of observ- 
ing that combination: 
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(17) 

Note that this sum naturally assigns more weight to more 
confident measurements. 

Assuming the lens and source redshift uncertainties 
can be added in quadrature, 

6E% = 5E% L +SS% S , (18) 

we find SEz ~ 0.177. 

Of course, these are just estimates for large ensembles. 
In practice, redshift probability distributions P(z) for 
individual galaxies will be properly folded into the P(T~c) 
determinations. Biased redshifts would yield biased Tc , 
the effects of which we study in Paper III. 

3.4. Projected Constraints from Large Surveys 

We now calculate the total uncertainty STc expected 
for large surveys with photometric redshifts. The com- 
bined lens model and time delay uncertainties are SEc T ~ 
0.344, based on extrapolation of the current empirical 
IQguril (|2007f ) finding ( %TTj) . We estimate uncertain- 
ties of SEz ~ 0.177 due to redshift uncertainties of 
Az L ~ 0.04(1 + z L ) and Az s ~ 0.10(1 + z s ) for the 
lenses and sources, respectively f ^3.3p . 

The simplest estimate of the total uncertainty is to 
add these uncertainties in quadrature: 

SE 2 = SE 2 Ct +SE%. (19) 

This yields 5E ~ 0.387. 

To be more precise, all of the uncertainties should be 
added in quadrature for each lens individually before 
combining them according to Eq. 1171 Repeating the anal- 
ysis in this way, we find SE ~ 0.402. 

Thus we expect large surveys with photometric uncer- 
tainties given above to yield STc ~ 40% /y/~N. We project 
ST C ~ 1-3% for PS1 (1,000 lenses) and ST C ~ 0.64% for 
LSST (4,000 lenses). 

Table[T]summarizes the progress we can expect to make 
in "Stages" corresponding t o those defined by the Dar k 
Energy Task Force (DETF; lAlbrecht et ail [20061 l2009f ): 
"Stage I" = current, "II" = ongoing, "III" = currently 
proposed, "IV" = large new mission. Again, we stress 
these are estimates of statistical uncertainties only. Large 
surveys are compared to dedicated monitoring and de- 
tailed analysis of a smaller sample of lenses. 

We might have made our analysis more sophisticated 
still, calculating SEc T , 5Ez L , and SEz s individually for 
each lens-source combination in our ensemble. Lenses 
and sources at higher redshift, for example, will be 
brighter and higher magnification cases on average, al- 
tering their 5Ec T somewhat. The approximations made 
in our above analysis should suffice for our purposes here. 

3.5. Quality vs. Quantity 

Thus far we have assumed that detailed observations 
and analysis would not be performed on the lenses. The 
alternative is to study fewer lenses in more detail, reduc- 
ing the uncertainties for each lens. In practice, we expect 



both strategies to be pursued and the combined power of 
both analyses to place the tightest possible constraints 
o n Tc- 

Moustakas ct al. (2008) have designed a mission con- 
cept that would be dedicated to monitoring a sample of 
four-image lenses, with the primary goal of constrain- 
ing fundamental properties of dark matter. This space- 
based Observatory for Multi-Epoch Gravitational Lens 
Astrophysics (OMEGA) would monitor 100 time delay 
lenses to achieve precise and accurate < 0.1 day time 
delay measurements. Supporting measurements would 
aim to reduce the model uncertainty of each lens to 5% 
(SEc = 0.05) and thus constrain Tc t o 5% with each lens , 
as claimed recently for B1608+656 (|Suvu et al.ll2009af) . 
These supporting measurements, including velocity dis- 
persion in the lens and characterization of the group envi- 
ronment (see discussion in £|6.2[) , would be carried out ei- 
ther with OMEGA itself or though coordinated efforts by 
ground-based telescopes and JWST. Spectroscopic red- 
shifts would also be obtained for the 100 lens galaxies 
and lensed quasars. 

Lenses targeted by OMEGA will be quads, enabling 
measurements of time delay ratios among the image 
pairs. This would provide co nstraints on the dark matter 
substructure mass function (|Keeton fc M oustakasl 120091 
iKeetonl 12009 , Moustakas et al., in preparation). 

Given lens models accurate to 5% for 100 galaxies, we 
might expect OMEGA to yield 5T C ~ 5%/VT00 = 0.5%. 
The time delays would be measured with sufficient pre- 
cision so as not to contribute significantly to the total 
uncertainty in STc- The multiple time delay measure- 
ments per lens (quad) also help reduce this contribution. 
Based on the expected time delay distribution for a sam- 
ple of quads (Paper I), we estimate that A(Ar) = 0.1- 
day uncertainties would inflate the Tc uncertainty only 
to - 0.515%. 

If both LSST and OMEGA obtain their measure- 
ments of Tc free of significant systematics, their com- 
bined power could further reduce the uncertainty to 
ST C ~ 0.4%. 

4. DEPENDENCE OF T c ON COSMOLOGY 

We expect LSST time delay lenses to constrain Tc to 
~ 0.64%. In this section we begin to explore how this 
"Stage IV" constraint translates to constraints on cos- 
mological parameters. We study the dependence of Tc 
on (h, fl m , fide, Ofc, wo, w Q ) for several cosmologies as out- 
lined in Table H 

4.1. Flat universe with a cosmological constant (h, 
f^A = 1 — ) 

First, we add a single free parameter Qa (in addition 
to h) in considering a flat universe with a cosmological 
constant (w = —1). Given STc = 0.64% from an ensem- 
ble with all lenses at zl =0.5 and all sources at zs = 2.0, 
we would obtain confidence contours shown in Fig. 2] 

The shape of these curves shifts somewhat as a function 
of zl and zs- Given an ensemble of lenses and sources 
with Gaussian redshift distributions zl = 0.5 ±0.15 and 
zs = 2.0 ±0.75 as discussed above, we begin to break the 
(h, Qa) degeneracy (Table EJ). Assuming a flat universe, 
Stage IV time delays could provide independent evidence 
for Da > 0. Whether this remains interesting by Stage 
IV remains to be seen. The constraints on h are certainly 
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Fig. 2. — Photometric redshift uncertainties' contributions to cosmological uncertainties in Tq. Left: Uncertainty in Tc (grayscale and 
contours) from lens redshift uncertainties of 0.04(1 + Zl,), plotted as a function of lens redshift zl and the lens-source redshift difference 
zg — zl. The dashed contours show the redshift distribution (1- and 2-er contours) assumed in this work. A dot at (z/,, Zs) = (0.5,2.0) 
marks the center of the distributions. Right: Same for source redshift uncertainties of 0.10(1 + z^). Note that the plots have different 
grayscales. For sources close to the lens (small zg — Zl), redshift uncertainties become catastrophic yielding large STc- Lens redshift 
uncertainties are also problematic at low Zj,. 



TABLE 1 

Estimated Current and Future Constraints on Tc 



Stage Experiment JVl quads Az a A(Ar) SE.c STc 



I 


current 


16 


38% 


spec 


2 days 


■ ■ ■ 8.6% 


11 


Pan-STARRS 1 


1,000 


19% 


phot 


2 days 


• • ■ 1.27% 


IV 


LSST 


4,000 


14% 


phot 


2 days 


0.64% 


IV 


OMEGA 


100 


100% 


spec 


0.1 day 


5% 0.5% 


IV 


LSST + OMEGA 










• • ■ 0.4% 



a Spectroscopic or photometric redshift measurements. For the latter we as- 
sume Az L = 0.04(1 + z L ) and Az s = 0.10(1 + z s ). 



TABLE 2 

Cosmologies explored in this work 



Cosmology 


h 








fide 1 ^A a 


Sly. 


wo 


w a 


Sect 


ions b 


Figures 


Flat universe with cosmological constant 


Free 


1 - 


-fi A 




Free (fi A ) 





-1 





i 


4.1 




HE 




Curved universe with cosmological constant 


Free 


1 - 


-(Ha 4 




Free (SI a) 


Free 


-1 







4.2 




IE 




Flat universe with constant w c 


Free 


1 - 


- f^de 


Free 





Free 







4.3 


El 


m 


EE2 


Flat universe with time- variable w 


Free 


1 - 






Free 





Free 


Free 


S 


4.4 


w 




General (curved with time-variable w) d 


Free 


1 - 






Free 


Free 


Free 


Free 


i 


5.3 




[SI [HE] 



Note. — We consider six cosmological parameters of which five are independent since Sl m + Sl^e + = 1- 
a When wo = — 1 and w a = 0, Side = ^Ai the cosmological constant. 

b In [0]the Tc dependencies are explored. In ij5] additional priors are assumed and time delay constraints are compared to those from other 
methods. 

c Given this cosmology, we assume a Planck prior in i|5.2l 

d Given a general cosmology, in i|5.3l we assume a prior of Planck + "Stage II" WL+SN+CL (see that section for details). 



tighter and would be improved by the introduction of a 
prior on f2 A , which we defer until fJSJ 

4.2. Curved universe with cosmological constant 
(h,Cl m ,ClA,Qk) 

If we relax the flatness parameter, adding another free 
parameter Q m (where curvature is determined by flk = 
1 — (f2 TO + n A )), we run into the degeneracy in Fig. [6l 
Plotted as colored bands are the (fi m , f2 A ) confidence 
contours assuming constant h — 0.7 given STc = 0.64% 
from an ensemble with all lenses and sources at zl, z$ = 
(0.5, 2.0). As h varies, these contours move as shown. 

An ensemble of lenses with a range of redshifts shrinks 



the confidence contours somewhat, as we see in Fig. [3 
though the strong (h, f2 m , fi A ) degeneracy remains. Even 
adopting an aggressive 3% prior on h, wc find neither f2 m 
nor J7 A can be constrained individually. However, the 
degeneracy does exhibit a strong preference toward a flat 
or nearly flat universe. Finally, we note the (h, fi m ,f2 A ) 
degeneracy can be more cleanly broken if our ensemble 
includes a significant fraction of lenses at zl = 1 and 
higher. 



4.3. Flat universe with constant dark energy EOS 
(h,tt de = 1-Ct m ,w) 



°- 4 rr 




1 10 100 1000 10000 

# lenses 



Fig. 3. — Constraints on STc as a function of ensemble size and 
observational uncertainties. The current ensemble has time delays 
measured to roughly A (At) = 2 day precision and spectroscopic 
redshifts measured for all lenses and sources. Future large surveys 
("quantity") should have similar time delay precisions but pho- 
tometric redshifts measured for lenses (Azr, = 0.04(1 + zr,)) and 
sources (A25 = 0.10(1 + 2g)). A dedicated campaign ("quality") 
could in principle obtain tight lens model constraints (<5H£ = 5%) 
with high-precision time delays (A(Ar) = 0.1 day) and spectro- 
scopic redshifts. 




h 



Fig. 4. — Confidence contours (1- and 2-cr colored bands) for 
(h, Q A = 1 - U m ) given "Stage IV" ST C = 0.64% obtained from 
an ensemble with all lenses and sources at z^, zg = (0.5, 2.0). 
Here we assume a flat universe with a cosmological constant (w = 
— 1). Also plotted are contours of constant 5c = Tc/Tc, true: where 
^C.truc ~ 0.99 for the input redshifts and cosmology. The input 
cosmology [h, f2 m , C2a) = (0.7, 0.3, 0.7) is marked with dotted lines 
and a white dot. 

Current cosmological constraints are consistent with a 
flat universe with a cosmological constant (as explored 
in §4.1|) . As a first perturbation to this model, it is com- 
mon to explore constraints onw^-1 while maintaining 
constant w in a flat universe. This cosmology has three 
free parameters (h, fi^cw) with f2 m = f — fide- 

Given enough data and appropriate priors, time delay 
lenses could place strong constraints on the dark energy 
equation of state parameter w (see fc|5. 2p . Figs. [5] and [5] 
explore the dependence of Tc on (w, fide) assuming a flat 
universe and constant w. 

4.4. Flat universe with time-variable dark energy EOS 
(h,fl de = 1 - fi m ,w ,w a/ ) 




h 



Fig. 5. — Confidence contours (1- and 2-cr colored bands) for (h, 
£2a = 1 — f2m) given STc = 0.64% and assuming a flat universe 
with a cosmological constant (w = —1). Each of the three fainter 
curves corresponds to all lenses and sources at the same pair of 
redshifts: z L , z s = (0.65, 2.75), (0.5, 2.0), (0.35, 1.25), as marked. 
Next we consider an ensemble of lenses and sources with Gaussian 
redshift distributions: z L , z s = (0.5±0.15, 2.0±0.75). These yield 
the tighter constraints (marked "ensemble"). The input cosmology 
(h,Q m ,QA) = (0.7, 0.3, 0.7) is marked with a white dot. 




Fig. 6. — Confidence contours (1- and 2-cr colored bands) for 
(f2 m ,OA) given STc — 0.64% obtained from an ensemble with all 
lenses and sources at z^, zg = (0.5, 2.0). The colored bands shift 
in (Q m ,f2A) space as h varies. A cosmological constant (w = — 1) 
is assumed. The input cosmology (h, Q m ,il\) = (0.7, 0.3, 0.7) is 
marked with a white dot. Flat cosmologies lie along the dotted 
line, and this line's intersection with the colored bands explains 
the strange shape of the colored bands in the previous plot. 

The most interesting constraints we can hope to place 
on dark energy are to verify or falsify the following: w = 
— f (cosmological constant) and w a = (constant w). In 
Fig. [10] we explore the dependence of Tc on (w , w a ) (see 
Eq. [TO]) . The colored bands are the constraints we could 
obtain given perfect knowledge of (h, fi m , fide)- The solid 
lines on the left show the curves' migration as a function 
of h. On the right, we also explore dependence on fide 
for a flat universe (il m + fide = 1)- 

5. COSMOLOGICAL CONSTRAINTS FROM FUTURE 
EXPERIMENTS 

We now consider the full parameter space 
(h, Q, m , fldei ^fcj w o, w a) and derive the constraints 
that may be placed on these parameters given con- 
straints on Tq along with various priors. Stage IV time 
delay constraints are compared to those expected from 
other experiments as estimated by the Dark Energy 
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(J) ensemble, h = 0.7 
fixed (zl, zs), h = 0.7 
ensemble, h = 0.7 ±—0 



0.6 



0.8 



1.0 



Fig. 7. — Additional confidence contours for (Q m ,Q\). The 
middle set of contours was plotted in the previous figure. The top 
set of contours assumes an ensemble of lenses and sources Zl,, zs 
= (0.5 ± 0.15, 2.0 ± 0.75). Finally, the bottom set of contours is 
for the ensemble and allowing a 3% uncertainty in h. 




Fig. 8. — Confidence contours (1- and 2-cr colored bands) for 
(w, Qd e = 1 — Q m ) assuming a flat universe with constant w given 
STc = 0.64% obtained from an ensemble with all lenses and sources 
at zl, zg = (0.5, 2.0). The colored bands shift in (w, fl^e) space 
as h varies. The input cosmology (/i, £7^, w) = (0.7, 0.7, -1) is 
marked with a white dot. 




0.0 



(•) ensemble, h = 0.7 
fixed (zl. zs). h m 0.7 
ensemble, h = 0.7 ±— £7 



1.0 



Fig. 9. — Confidence contours for (w, Q^e = 1 — f2 m ), assuming 
a flat universe. As in Fig. [7] we plot a "Stage IV" ensemble of 
lenses at a range of redshifts, the lenses all at the same redshift, 
and the ensemble allowing 3% uncertainty in h. 



TABLE 3 

Stage IV Fisher matrix expectation for (h, 0.^,0,^, wo,w a ) 
given 8T C = 0.64% 





h 




Q k 


Wo 


w a 


h 


49824.9224 


-1829.7018 


-4434.2995 


4546.8899 


122.5319 




-1829.7018 


88.3760 


200.9795 


-189.2658 


-8.4386 




-4434.2995 


200.9795 


463.5732 


-445.5690 


-17.9694 


WO 


4546.8899 


-189.2658 


-445.5690 


441.9725 


15.2981 


Wo 


122.5319 


-8.4386 


-17.9694 


15.2981 


1.0394 



Task Force (|Albrecht et all l2fM l2009t ). To efficiently 
explore this parameter space, we perform Fisher matrix 
analyses. 

5.1. Fisher Matrix Analysis 

The Fisher matrix formalism provides a simple way to 
study uncertainties of many correlated parameters. Con- 
straints from various experiments and/or specific priors 
may be combined with ease. A "quick-start" instruc- 
tional guide and software are provided in a companion 
paper (jCoel 120091 ). Fisher matrices approximate all un- 
certainties as Gaussians. The true uncertainties may be 
somewhat higher and non-Gaussian. The full informa- 
tion of the dependencies as shown in [j4] is not retained. 
Yet as cosmological parameters are constrained close to 
their true values, these approximations should suffice. 

As above we consider a "Stage IV" ensemble of time 
delays which constrains Tq to 0.64% with Gaussian dis- 
tributions of lens and source redshifts (zl — 0.5 ± 0.15; 
zs = 2.0±0.75). Assuming such a Gaussian distribution 
for 7c and the aforementioned redshift ensemble, we cal- 
culate (numerically) the Fisher matrix for cosmological 
parameters of interest. The Fisher matrix consists of 
partial derivatives of x 2 with respect to the parameters. 
For parameters (j>i,Pj), element in the Fisher 

matrix is given by 



1 d X 2 



2 dpidpj 



(20) 



The Stage IV (ST C = 0.64%) time delay Fisher ma- 
trix is given in Table [3] for the cosmological parameters 
(/i, fide, fifc, wo, w a ). The Fisher matrix may be easily 
scaled to other 67~c values. For example, to scale from 
LSST (4,000 lenses; 5T C = 0.64%) to Pan-STARRS 1 
(1,000 lenses; 8Tq — 1.27%), simply divide all the values 
in the Fisher matrix by 4 = 4, 000/1, 000 = (1.27/0.64) 2 . 
Or multiply them by 1.6 = (0.64/0.4) 2 to explore the 
LSST + OMEGA constraints (6T e = 0.4%). If one 
is interested in constraints on O m = 1 — (Qde + ^/c), 
uj m = Q m h 2 , or any other related variable, a tran sfor- 
matio n of variables can be performed as outlined in iCod 
((20091) . 

In Fig.[Tl]we show the time delay constraints possible 
on all parameters and pairs of parameters assuming per- 
fect knowledge of all the other parameters. These plots 
can be compared to those presented in fJU Such perfect 
priors are unrealistic, but they help to demonstrate the 
parameter dependencies and degeneracies. 

5.2. Flat universe with constant w 

We first consider the simple case of a flat universe with 
constant w. This is a common perturbation to the con- 
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Fig. 10. — Left: Confidence contours (1- and 2-a colored bands) for (wo,w n ) given STg = 0.64% obtained from an ensemble with all 
lenses and sources at Zl> zg = (0.5, 2.0) and assuming h = 0.7 and perfect knowledge of (f2 m , f!d e ). As shown, these bands shift in (wo, w a ) 
space as h varies. The input cosmology (/i,wo,w ) = (0.7, -1, 0) is marked with a white dot. Right: Dependence of the (wo,w a ) contours 
on O^e, assuming a flat cosmology. Dashed lines show the h dependence from the left plot. Solid lines of increasing thickness show contours 
of Sljj,, decreasing in 0.1 increments. 




h n o 



Fig. 11. — Constraints placed on pairs of parameters derived from our Fisher matrix analysis assuming perfect knowledge of all other 
parameters given &Tc = 0.64% obtained from an ensemble with Gaussian distributions of lens and source redshifts (zl = 0.5 ± 0.15; 
zg = 2.0 ± 0.75). All off-diagonal plots show 1- and 2-cr colored ellipses. Along the diagonal are constraints on individual parameters 
assuming perfect knowledge of all others. The y axes along the diagonal are units of relative probability, different from the off-diagonal 
plots. 
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cordance cosmology. The goal is to detect deviation from 
w = — 1, equivalent to the cosmological constant A. This 
3-parameter cosmology (h, fide-, w, with fl m = 1 — fide) 
was explored above in ^4.31 

The top row of Fig. [T2] shows Stage IV time delay con- 
straints with a Planck prior in a flat universe with con- 
stant w. Given these priors, we estimate that time delays 
will constrain h to 0.007 (~ 1%), fide to 0.005, and w to 
0.026 (all l-cr precisions). 

In the bottom row of Fig. Q21 we compare these time 
delay constraints (TD) to those expected from other 
methods: weak lensing (WL), baryon acoustic oscilla- 
tions (BAO), supernovae (SN), and cluster counts (CL). 
We consider "optimistic Stage IV" expectations from 
these methods as calculated by the Dark E nergy Task 
Force (DETF; lAlbrecht etHI I2006L I2009T ) and made 
available in the software DETFast 9 . A Planck prior (also 
calculated by the DETF) is again assumed for all exper- 
iments. 

In manipulating the DETF Fisher matrices we adopt 
their cosmology (fl m , fide, h) = (0.27, 0.73, 0.72), but we 
revert to our chosen cosmology (fl m , fide, h) = (0.3, 0.7, 
0.7) for the rest of our analysis. These differences have 
n egligible impact on ou r results. 

iLewis fc Ibatal (|2002|1 considered similar constraints 
from time delay lenses but found much weaker con- 
straints on (h, w), even with all other cosmological pa- 
rameters fixed. One of the cases they considered was 
500 lenses with 15% uncertainty each, which translates 
to 15%/V500 = 0.66% total uncertainty, very similar to 
the 0.64% uncertainty we estimate for LSST given 4,000 
lenses with a much higher uncertainty (effectively 40%) 
assumed per lens. For this case, they find 0.99 < h < 
1.10 and -1.48 < w < -0.88 (95% confidence). When 
we perform a similar analysis, assuming STc — 0.64% 
and perfect knowledge of (f2 m , fide, ^k, w a ), we obtain 
similar uncertainties (without biases, by construction): 
h = 0.7 ± 0.02 and w = -1 ± 0.21 (l-cr). But with the 
addition of a Planck prior, even while relaxing the perfect 
prior on (fl m , fide, flk,~Wa), we find improved constraints 
of h = 0.7 ± 0.007 and w = -1 ± 0.026 (l-cr). Planck 
clearly complements the strong lensing constraints well 
to produce tight constraints on (h, w). 

5.3. General Cosmology 

We now assume a general cosmology allowing for cur- 
vature and a time-varying w. To help constrain this 
larger parameter space (h, fide, flk, wo, w a , with fl m = 
1 — (fide + flk)), we add additional priors. In addition to 
the Planck prior, we adopt "Stage II" (near-future) con- 
straints from weak lensing (WL) + supernovae (SN) + 
cluster counts (CL), all as calculated by the DETF. The 
DETF uses this prior (in addition to Planck) in many 
of their calculations comparing the performance of Stage 
III - IV techniques. 

The Stage II DETF WL + SN + CL prior yields the fol- 
lowing uncertainties: Ah = 0.031 (4.4%), Afl de = 0.023, 
Afl k = 0.010, Aw = 0.128, Aw a = 0.767 (along 
with various covariances between parameters). The ad- 
dition of the Planck prior reduces these to: Ah = 0.017 
(2.4%), Aflde = 0.012, Afl k = 0.003, Aw = 0.115, 

9 http: //www .physics .ucdavis . edu/DETFast/ 



Aw a = 0.525. Note that Stage II WL+SN+CL con- 
strains h well enough (to 4.4%) that an HST Key Project 
prior (h — 0.72 ± 0.08) appears to be unnecessary. Even 
SHOES (h = 0.742 ± 0.036, or 4.9%) provides a weaker 
constraint on h. However, as noted in the introduction, 
these combined WL+SN+CL experiments yield a predic- 
tion of h based on an assumed cosmological model and 
are n o substitute for local measurements of h ()Riess et al.l 
f2009h . 

These Stage II constraints are also rather optimistically 
combined, assuming that all experiments have converged 
on the same best fit cosmology without systematic offsets 
among them. The true Stage II constraints should be 
somewhat weaker. 

Plotted in Fig. 1131 are time delay constraints assuming 
a prior of Planck + Stage II WL+SN+CL. A progression 
is shown from Stage I (present) time delay constraints 
(ST C = 8.6%) through Stage II (6T C = 1.27%) and on to 
Stage IV (STc = 0.64%). The current constraints barely 
improve upon this aggressive prior. While the Stage II 
- IV constraints certainly improve upon the prior, note 
that the outer bounds of the time delay and prior ellipses 
nearly intersect. This indicates that the size of the time 
delay ellipse is controlled by that of the prior, at least 
for these constraints and prior. Were the prior signifi- 
cantly weaker or the time delay constraints significantly 
stronger, we have verified that the time delay ellipses 
would shrink well within the prior ellipses. 

In Fig. [14] we compare Stage IV time delay constraints 
to those expected from other methods for various pa- 
rameters of interest. Plotted are constraints on (h, flk), 
(h, wq), and (wo, flk), and (wo,w ). An example of how 
these constraints combine is given in £15.4.21 

We give extra attention to constraints on the dark en- 
ergy parameters (wo,w a ). The DETF figure of merit 
(FOM) for a given experiment is defined as the inverse 
of the area of the ellipse in the (wo, w a ) plane. In Fig. [T5] 
we plot FOM for various experiments versus the "pivot 
redshift", defined as follows. For a time-varying w(z), 
time delays constrain w best at z s a 0.31. This red- 
shift is known as the pivot redshift (jHuterer fc Turner! 
l200H:IHu fc Jah3[200l and can also be calculated sim- 
ply from the (w ,w a ) constraints (Coc 2009). As in the 
previous plot, we assume a prior of Planck + Stage II 
(WL+SN+CL). 

5.4. Time delays do not simply constrain h 
5.4.1. Relaxing the "perfect prior" on (flm, fide, flk, wo,w a ) 

To date, analyses of time delay lenses have quoted un- 
certainties on Tc as uncertainties on h, assuming Sh = 
STc- This assumption has been valid to date, but future 
constraints on h will be weaker than the constraints on 
7c, that is Sh > ST C . 

This is demonstrated in Fig. [TB] left. The dashed 
line shows Sh = STc, or the "perfect prior" on 
(fl m , fide, ^k, wo, w a ) generally assumed in analyses. For 
future samples (at the left side of the plot), as this 
prior is loosened, we find Sh > i57c- In Fig. [TB] right, 
we plot Sh/STc- For example, given a "Stage II" prior 
on WL+SN+CL, and LSST constraints on time delays 
(STc = 0.64%), we find Sh ~ 2.2ST C ~ 1.4%. Alterna- 
tively, assuming a Planck prior in a flat universe with 
constant w, we would find Sh ~ lASTc ~ 0.90%. 
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Fig. 12. — Top row: Cosmological constraints from "Stage IV" time delays plus a Planck prior in a flat universe with constant w. We 
assume an ensemble of time delays which constrains T~c to 0.64% (see text for details). Time delays plus Planck constrain h to 0.007 (1%), 
Qd e to 0.005, and w to 0.026 (all 1-cr precisions). Bottom row: Comparison of "optimistic Stage IV" constraints expected from time delays 
(TD), weak lensing (WL), supernovae (SN), baryon acoustic oscillations (BAO), and cluster counts (CL). The time delay constraints are 
as plotted in the top row. For the other experiments we use Fisher matrix calculations provided by the Dark Energy Task Force (DETF). 
For each parameter pair, experiments are plotted in order of FOM cx (Ellipse Area) -1 , with the best experiment on top. 



(Note that the Stage II WL+SN+CL prior claims a 
constraint of 5h = 0.03, such that it outperforms current 
constraints from time delays 5h = STq.) 

5.4.2. Time delays provide more than constraints on h 

In the introduction we commented on the ability of 
any experiment to improve constraints on w and £lk sim- 
ply by tightening the constraints on h. Several methods 
ha ye the poten tial to further improve the constraints on 
h (Oiling 2007). Do time delays offer more than a simple 
constraint on h for the purposes of constraining the dark 
energy equation of state? 

In Fig. [TT1 we compare Stage IV time delays (left) to 
a simple h constraint (right) in ability to constrain dark 
energy. Each is combined with Stage IV supernovae con- 
straints plus a prior of Planck + Stage II WL+SN+CL. 10 
We find time delays are more powerful than the simple 
h constraint. The (SN + TD + prior) figure of merit 
(FOM) on (w ,w Q ) is ~ 19% higher than that from (SN 
+ H + prior). 

The "H" constraint Sh = 0.009 was chosen such that 
when combined with the prior, the resulting Sh would 
equal that from TD + prior. Both H + prior and TD + 
prior yield Sh — 0.008. However we find TD outperforms 
even a perfect H prior (Sh ~ 0) by 13%. Simply put, 
the time delay constraints on (f2 m , Side, wo, w a ) are 
clearly making contributions. 

When combined with experiments other than SN, TD 
offers less marked improvements over H constraints. Re- 

10 Strictly speaking we have not taken the proper care in com- 
bining constraints from the Stage II and Stage IV supernova ex- 
periments, as their nuisance parameters have been marginalized 
over in the DETF Fisher matrices. But this analysis will suffice for 
illustrative purposes here. 



placing SN with BAO, WL, and CL, we find TD outper- 
forms H by 7%, 5%, and 3%, respectively. 

5.5. Lens and Source Redshift Distribution 

We have been considering the Gaussian redshift distri- 
butions zl = 0-5 ± 0.15, z$ = 2.0 ± 0.75 introduced by 
iDobke et al.l ([2009:) as reasonable approximate assump- 
tions for near-future missions. We find that the cosmo- 
logical parameter constraints are not extremely sensitive 
to variations in these redshift distributions. 

For 5T C = 0.64% plus our Planck + Stage II 
(WL+SN+CL) prior, we find the following. A lower 
tighter lens redshift distribution of Zl = 0.2 ±0.1 im- 
proves the constraint on h by 22% and on Vide by 12% 
at the expense of the wo and w constraints, which de- 
grade by 8% and 10%, respectively. A higher tighter lens 
redshift distribution of zl — 1.0 ± 0.1 has less leverage, 
as the h and f2<j e constraints degrade by 15% and 14%, 
respectively with little benefit to the other parameters. 
Neither broader lens redshift distributions nor variations 
on the source redshift distribution have much impact on 
the parameter constraints. 

When time delay constraints are tighter (STc < 
0.64%), with the same priors, the lens redshift distribu- 
tion begins to have a greater impact. We reserve study 
of such "beyond Stage IV" constraints for future work. 

6. SYSTEM ATICS 

As with any measurement, there are many potential 
sources of systematic bias, as alluded to throughout this 
work. At the risk of putting the cart before the horse, we 
have presented systematic-free projections for time delay 
cosmological constraints. These should serve to motivate 
a more considered look at systematics, in the context of 
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Fig. 13. — Cosmological constraints from time delays in a general cosmology assuming priors of Planck + Stage II (WL+SN+CL) 
as calculated by the DETF. A progression is shown from the prior (outermost ellipse, 2-er) to Stage I (current) time delay constraints 
(STc = 8.6%; gray ellipse, 2-cr) to Stage II constraints (STc = 1.4%; black ellipse, 2-<r) to Stage IV constraints (STq = 0.64%; colored 
ellipses, 1-cr and 2-cr). Along the diagonal are plotted constraints on individual parameters marginalizing over all others. The y axes along 
the diagon al a re units of relative probability, different from the off-diagonal plots. For each plot, we marginalize over all other parameters, 
unlike Fig. 1111 in which unrealistic perfect priors were assumed for illustrative purposes. 



the behavior of random uncertainties in these studies. 
Ideally, efforts should be undertaken to reduce system- 
atics on a timescale comparable to that presented here 
(e.g., 0.64% by "Stage IV"). If this cannot be accom- 
plished, we study prospects for estimating cosmological 
parameter s in spite of large residual systematic biases in 
Paper III (jCoe fc Moustakasll2009cl ). 

Here we discuss briefly the greatest potential sources of 
systematic bias. We should consider which of our main 
sources of statistical uncertainty (lens modelling, redshift 
measurements, and time delay measurements) could also 
contribute significant systematic bias. Time delay un- 
certainties are generally not expected to be biased in 
any preferred direction. Redshift biases are somewhat 
worrisome but will not be discussed further here. Most 
daunting are potential biases due to imprecise lens mod- 
eling. 

Whether we determine the appropriate lens model for 
the "typical" ( "average" ) lens in an ensemble or we con- 
strain each individual lens model well, we must use the 
following tools to measure lens properties. The largest 
statistical uncertainties and potential systematic biases 
involve measurements of the lens mass density slope and 
perturbing mass sheets. 



6.1. Lens Mass Density Slope 



Regarding mass slope, this paper has focused on the 
statistical strategy which assumes that we know the cor- 
rect mean of mass slopes. Evidence currently suggests 
that lenses are isothermal (a = 1, 7 = 2) 11 on average. 
Yet a recent analysis of 58 SLACS lenses finds a slightly 
higher average slope of 7 = 2.085±o;oig(stat.)±0.1(syst.) 
(|Koopmans et al. 2009b). If the average proved to be ex- 
actly 7 = 2.085, this would result in an 8.5% bias in Tq 
(STc = S"//2 = 8a) were we to assume an average of 
7 = 2 instead. 

Mass profile slopes for individual lenses are determined 
by measuring mass within two radii: the Einstein radius 
(from the positions of multiple images) and a smaller 
radius (from velocity dis persions). The latter req uire de- 
tailed spectroscopy (e.g.. iKoopmans et al.ll2006h . It will 
not be feasible to obtain the required measurements for 
all time delay lenses detected in future surveys, but small 
samples of these can be selected for such detailed study. 

6.2. Mass Sheets 

Mass sheets can be equally harmful as a source of sys- 
tematics as 7c bias also scales linearly with projected 

We use two definitions common in the literature regarding 
lens slope: two-dimensional mass surface density ft oc r 2-Q , and 
three-dimensional mass surface density p o ar These parameter s 
are related by a + 7 ss 3 (see discussion in van de Veil et al. 2009). 
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Fig. 14. — Comparisons of "Stage IV" constraints possible from time delays (TD), weak lensing (WL), supernovae (SN), baryon acoustic 
oscillations (BAO), and cluster counts (CL) in a general cosmology (allowing for curvature and a time-variable w). For TD, we assume an 
ensemble which constrains T~c to 0.64% (see text for details). For the rest we use "optimistic Stage IV" expectations calculated from Fisher 
matrices provided by the Dark Energy Task Force (DETF). A prior of Planck + Stage II (WL+SN+CL) is assumed for all five experiments 
and is plotted in gray. For each parameter pair, experiments are plotted in order of FOM oc (Ellipse Area) -1 , with the best experiment on 
top. 
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Fig. 15. — Dark energy figure of merit (FOM oc 
((wfj,w a ) Ellipse Area) -1 , normalized relative to the prior) versus 
pivot redshift for various "optimistic Stage IV" experiments with 
a prior of Planck + Stage II (WL+SN+CL) The pivot redshift is 
the redshift at which w(z) is best constrained. 

mass density, 5Tq ~ k. Mass sheets can result from 
both mass within the lens group environment and mass 
along the line of sight (over- or under-densities) all the 
way from source to obse rver. The former is th e domi- 
nant effect. Simulations (|Dalal fe Watsonl 120051 ) suggest 



that group members contribute K C nv = 0.03 ± 0.6 dex 
(i.e., log 10 (K cnv ) = log 10 (0.03) ± 0.6) for a l-a upper 
bound of K cnv = 0.12, or 12% bias on Tq. Mass along 
the line of sight is generally lower and more nearly fluc- 
tuates about the cosmic average but should also be ac- 
counted for. iHilbert et al.l (|2007l ) measured mass along 
the lines of sight to strong lenses in the Millennium sim- 
ulation. For sources at z$ = 2, the central 68% span 
-0.0355 < ki os < 0.0475 (Paper I). 

Efforts are made to measure K onv for individual 
lenses via spectroscopic (and photom etric) studies (e.g., 
iMomcheva et al.l 120061: lAugerl 120081 ) and simulations 
which estimate the effects of nearby neighbor s (e.g., 
IKeeton fc Zabludoffll200l IDalal fc Watsonll2005D . Simi- 
lar studies also attempt to identify groups along the line 
of sight and estimate t heir mass sheet contributions (e.g., 
iFassnacht et al.ll2006D . 

The alternative is a statistical approach. Measure- 
ments of K cnv or Ki os would not be required for individual 
lenses if we had knowledge of the distributions P{n env ) 
and P(ki os ) for strong lenses. These distributions could 
be obtained from simulations, and one could attempt 
to correct for the ex pected bias for lenses to reside in 
high density regions (Dal ai fe Watsonl 120051 : lOguri et all 
[20051) . However, one might wonder whether these dis- 
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Fig. 16. — Demonstration that Sh > STc for future ensembles. Left: Constraints on h versus constraints on T~c for various priors. Along 
the top horizontal axis we plot experiments with corresponding STc: current constraints (8.6%), Pan-STARRS 1 (1.27%), LSST (0.64%), 
OMEGA (0.5%), and LSST + OMEGA (0.4%). The priors are different combinations of the following: Planck, a flat universe, constant w, 
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Fig. 17. — Left: Combined constraints on (wq, w a ) from Stage IV time delays (TD) and supernovae (SN). A prior of Planck + Stage II 
(WL+SN+CL) is assumed. The TD + prior constraint yields Sh = 0.008 (not shown). Right: Similar plot combining Stage IV SN with a 
Sh = 0.009 constraint on Hubble's constant (that which also yields Sh = 0.008 when combined with the prior). Time delays yield a 19% 
improvement in figure of merit (FOM oc ((wo,w„) Ellipse Area) -1 ), versus the constraint on h alone. SN + TD shows the most dramatic 
such improvement vs. SN + H. Replacing SN with the other experiments (BAO, WL, CL) we find lesser improvements vs. H of 7%, 5%, 
and 3%, respectively. 



tributions and corrections would prove accurate to the 
percent level. Any errors would yield residual systemat- 
ics in our estimation of Tc- 

To aid such a statistical approach, lenses in obvious 
groups can be excluded from the analysis leaving only 
those systems with low K env . Such low mass systems 
would introduce smaller biases, though a detailed explo- 
ration of this approach will await future work. 

7. CONCLUSIONS 

We have presented the first analysis of the potential of 
gravitational lens time delays to constrain a broad range 
of cosmological parameters. The cosmological constrain- 
ing power ST C was calculated for Pan-STARRS 1, LSST, 
and OMEGA based on expected numbers of lenses (in- 
cluding the quad-to-double ratio) as well as the expected 
uncertainties in lens models, photometric redshifts, and 
time delays. Our Fisher matrix results are provided to 
allow time delay constraints to be easily combined with 
and compared to constraints from other methods. 

We concentrate on "Stage IV" constraints from LSST. 



In a flat universe with constant w including a Planck 
prior, LSST time delay measurements for ~ 4, 000 lenses 
should constrain h to ~ 0.007 (~ 1%), n de to - 0.005, 
and w to ~ 0.026 (all 1-a precisions). We compare 
these results as well as those for a general cosmology to 
other "optimistic Stage IV" constraints expected from 
weak lensing, supernovae, baryon acoustic oscillations, 
and cluster counts, as calculated by the Dark Energy 
Task Force (DETF). 

Combined with appropriate priors (those adopted by 
the DETF) , time delays provide modest constraints on a 
time-varying w(z) that complement the constraints ex- 
pected from other methods. Time delays constrain w 
best at zkj 0.31, the "pivot redshift" for this method. 

We find that LSST and OMEGA represent about an 
even trade in "quantity versus quality" in terms of con- 
straining cosmology with time delays. LSST could yield 
(57c ~ 0.64% by measuring time delays for 4,000 lenses, 
while OMEGA could yield 6T C ~ 0.5% by obtaining 
high-precision time delay measurements and lens model 
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constraints for 100 lenses with spectroscopic redshifts. 
The combined statistical power of these two missions 
could further improve the cosmological constraints to 
5T C ~ 0.4%. 
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